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Abstract 

The energy range and quality of strong-interaction data from recent years demand the 
study of higher orders in perturbative QCD, and of nonperturbative effects. I discuss a 
selection of recent progress in the theory of QCD at high energy, including examples from 
perturbative resummation, nonperturbative power corrections and the Regge limit. In each 
case, techniques of factorization play a central role. 



1 Introduction 

Current studies in QCD are motivated partly by its importance in the production and detection of 
new physics, but also, and in very large part, by the challenges of quantum chromodynamics itself. 
QCD may be thought of as an exemplary quantum field theory, exhibiting asymptotic freedom, 
confinement, chiral symmetry breaking and so on, the workings of which are all available, and 
unavoidable, within present energy ranges. It is a vast subject, encompassing perturbative, heavy- 
quark, nonrelativistic and lattice QCD, all the way to nuclear physics. In a very real sense, what 
are sometimes called "tests of QCD" are tests of quantum field theory itself. 

1 Based on a talk presented at the IVth Rencontres du Vietnam, International Conference on Physics at Extreme 
Energies, Hanoi, Vietnam, July 19-25, 2000. 
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Figure 1: DIS cross section for neutral and charged currents at high momentum transfers. Q 

In this talk, I will discuss some recent efforts to narrow the gap between the high-energy, par- 
tonic and low-energy, hadronic descriptions of QCD, starting from the high-energy side, through 
the study of higher-order corrections in perturbation theory and of nonperturbative power correc- 
tions. I will highlight the use of QCD factorization as an organizing principle in these investiga- 
tions. 

The past decade has been a proverbial golden age of hadronic data at high energy, in terms 
both of coverage and of quality. There is no room here to do justice to the data revolution of 
the 1990's, the work of a generation of accelerators that has reached maturity: LEP, HERA and 
the Tevatron. Each has produced spectacular successes for our current picture of QCD, but each 
has provided its share of puzzles. Two examples must suffice. Our control of strong interaction 
corrections in inclusive cross sections such as deep-inelastic scattering (DIS) may be gauged from 
Fig. H, showing data from HERA. []J Here, the total ep cross sections at momentum transfers 
from tens into hundreds of GeV track the standard model predictions, which include extensive 
input from perturbative QCD evolution, described below. 

On the other hand, Fig. || shows the transverse momentum distribution for b-quark production 
at the Tevatron. |2| This data is typical of cases where our present theory is partly adequate, 
partly not. Overall, the theoretical px spectrum has the correct shape, but the normalization of 
the theory is too low at low pr, even as it gradually approaches the data at larger pt- Much of the 
effort described below is aimed at using such apparent discrepancies as guides to the perturbative 
and nonperturbative structure of the theory. To see how, let us briefly review the elements of our 
present theoretical framework for high-energy QCD, based on factorization. 

2 The Unity of QCD Factorizations 

The application of perturbation theory to high energy QCD begins with asymptotic freedom and 
infrared safety, illustrated by the perturbative expansions of the total cross sections for e + e~ 
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Figure 2: Transverse-momentum distribution for b-quark production at the Tevatron. 

annihilation to hadrons, and to final-state jets: 

Q 2 a(Q 2 , /i 2 , aM) = £ c n (Q 2 /fi 2 ) <(//) + 0(l/Q p ) = £ c n (l) a n s {p) + 0(1/ (?) , (1) 

n n 

where the c n are dimensionless coefficients. For an infrared safe cross section, the c n are free of 
dependence on fixed mass scales (such as light quark masses), which are absorbed into corrections 
that are suppressed by some power, p, of the c.m.s. energy, Q. In most cases, the c n are fully known 
only for a few, low orders. Because the cross section is a physical observable, it is independent 
of the renormalization scale, /i, which can therefore be chosen to equal Q. For an asymptotically 
free theory, the larger is Q, the better any finite-order approximation becomes. 

Relatively few cross sections are quite this simple, however, but whenever a reaction involves 
a scattering at large momentum transfer, or the creation or decay of a heavy state, we may 
isolate its short- distance components, which can be treated perturbatively, from its long-distance, 
nonperturbative components. This is a procedure known as factorization, which generalizes the 
operator product expansion. 

For a factorized cross section, Eq. (|l|) is replaced by an expression of the general form, 

Q 2 a(Q, x) = lu(Q/h, x/£, a,(fi)) ® /(£, n) + 0(l/Q?) , (2) 

with a " hard-scattering", or coefficient, function u, which is short-distance and perturbative, in 
convolution with a "soft" function /, which is long-distance and nonperturbative. In DIS, with q 
the momentum transfer, x = 2p-q/Q 2 , but more generally it represents any dimensionless ratios of 
large momentum scales. The dimensional variable \x is the factorization scale, separating long and 
short distances. As in (|l|), the physical cross section is independent of fi. In DIS, the soft function 
/ is a parton distribution function (PDF). For hadron-hadron scattering, we have two PDF's in 
convolution form. In these cases, the convolution in (0) is in terms of fractional momenta, £, of 
the (one or more) partons that initiate the hard-scattering process. Factorization is more general 
than this, however, and we shall encounter other examples below. Nearly always, the soft function 
can be interpreted as the matrix element of some (usually nonlocal) operator in QCD. 

The basis of factorization is always the quantum-mechanical incoherence of dynamics at very 
short distances from that at long distances, and, in Minkowski space, the mutual incoherence of 
the dynamics of particles whose relative velocity approaches the speed of light. 
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Whenever there is factorization, there is evolution, a consequence of the independence of the 
physical cross section from the factorization scale, 



l^-^lna(Q,x,m) = //^- In {u(Q/n, x/f, «,(//)) ® /(£,//)} = 0. (3) 

Because / and a> have in common only the parton momentum fraction and a s , separation-of- 
variable arguments imply complementary equations for / and to: 

\i — = P{z,a s ) <g> f{£/z, n) and lo[Q n,r)z,a s ) ® P{z,a s ) = -// , (4) 

a/i a/i 

in terms of convolutions with splitting functions P(z,a s ). The first of these "DGLAP evolution" 
relations enables us to take PDFs determined at some reference scale, and extrapolate to 
higher, or lower, scales, wherever the running coupling is not too large. 

Measurements of the strong coupling based on these methods give a s [Mz) ~ 0.12, which 
suggests that at around 100 GeV, 0(a 2 s ) is about one percent. This is the nominal level of accuracy 
to which perturbative QCD may aspire at "next-to-next-to-leading" (NNLO) order. Referring to 
Eq. (HD, the c n 's are known to NLO (n = 1) for cross sections with up to four jets in e + e". [||] 
To date, complete NNLO calculations are available only for one-scale problems: the total e + e~ 
cross section, DIS and Drell-Yan. Two loops are the current frontier for finite-order perturbative 
QCD, and the past year has seen significant progress toward the exact computation of two-loop 
scattering amplitudes and coefficient functions to. Q At the same time, to use two-loop coefficient 
functions, it will be necessary to have the splitting functions P at three loops; and here also 
important progress has been reported within just the past few months. || [| 

Beyond exact calculations, DGLAP evolution is the first among a set of methods that enable 
us to probe properties of QCD perturbation theory to all orders. Each of these methods is based 
upon the separation of dynamics at different length scales. 



3 Resummation for Inclusive and Exclusive Cross Sections 

Hard-scattering, such as jet, heavy quark and high-p T photon production, depends on a complex 
combination of the long-distance dynamics of the external hadrons, the short-distance perturbative 
subprocess and the properties of final states included in the cross section. Using concepts of 
factorization, we are learning to treat a widening variety of such cross sections, and increasingly 
to compute classes of higher-order corrections. We describe below three applications of current 
interest. 



3.1 Partonic Threshold 

Our first example is so-called threshold resummation, which applies to inclusive hard-scattering 
hadronic cross sections AB — * F+X, with F=7*, W, Z, jets, heavy quark, etc. of invariant mass Q. 
We are interested in higher-order corrections to the perturbative hard-scattering function uo a ^F 
in factorized cross sections, Eq. fl2|), at partonic threshold, z = Q 2 /s — > 1, where the partons a and 
b have just enough c.m.s. energy to produce the observed final state. In the threshold region, 
there is an incomplete cancellation of emission and virtual radiative corrections, which leads to 
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singular corrections of the form 



(r) 



7r / r! 



ln 2r - 1 (l 



(5) 



at rth order, with C q = Cp, C g = Ca- In effect, for z — > 1, there are two hard scales, Q and 
(1 — z)Q. Because these singular distributions arise from soft-gluon radiation, sensitive to the 
scale (1 — z)Q, but not to the scale Q, they may be "refactorized" from the hard-scattering into 
functions that depend only on the flow of color at a truly short- distance hard-scattering, which is 
sensitive only to Q. As for Eq. (||), the new factorization implies a new evolution equation, from 
which we derive a resummation in moment space, 



U)ab^F(Q/^,N) 



1 

exp 



dzz N 1 uj ab ^ F (Q//j,,z) 



1 Z N-1 

dz 

o 1 — 



1 



dm 2 



z J(i-z) 2 q 2 m z 



A ab (a s (m)) 



(6) 



where A ab = {C a + C b )(a s /ii) + ... is an expansion in a s . 

Threshold resummation is currently being explored for most of the basic inclusive cross sections. 

|8| An important consequence is a reduction of factorization-scale dependence in resummed 
cross sections. In NLO calculations, dependence on /i begins at the next order, 0(a^), but often 
this residual sensitivity is uncomfortably large. In the resummed function Eq. (||), however, the 
/i-dependence is determined by 



din u ab (N,n) 
din ft 



A ab (a s ) InN 



d\n[f a/A (N,Li)f b/B (N,ii)} 
(iln/x 



(7) 



where to the right the /'s are moments of the parton distributions. Factorization requires that 
the function A ab that appears in Eq. (^|) is exactly the same as the sum of the IniV terms in 
the moments of the splitting functions P aa and P bb . This leads to a very significant decrease 
in sensitivity to the factorization scale compared to previous, fixed order calculations. [10, |TT 



This is only the beginning of applications of threshold resummation, however, and we anticipate 
important applications to the determination of parton distributions and to the improvement of 
predictions for new particle production. 



3.2 Power Corrections: Universality and Beyond 

Jet cross sections in e + e~ annihilation are defined by adjustable parameters, whose variation 
mediates between fully inclusive and nearly exclusive cross sections. As such, they are ideal for 
testing and improving our understanding of QCD at intermediate distances. 

The most-studied examples involve light-mass dijet pairs. Dijet events, which dominate the 
annihilation cross section at high energy, can be described in terms of event shapes. Perhaps the 
best known of these is the thrust, T. The thrust of an e + e~ event is determined approximately 
by finding an axis that maximizes the quantity T = 1 — (mf + m|)/Q 2 , where m x and m 2 are the 
invariant masses of the sums of all particle momenta within the two hemispheres defined by this 
axis. As T — >• 1, the final state is characterized by two well-collimated jets. A number of other 
familiar event shapes may be derived from jet masses in a similar way. 
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For light-mass dijets, the relative velocities of the two jets again insures that their dynamical 
developments into the final state are mutually independent, This independence results in factor- 
ization at the level of cross sections, and leads, in much the same manner as above, to evolution 
equations, and to T-moments that are quite similar to the z-moments of threshold resummation 
in Eq. (1), 



/ 



dTT 



N-l 



rfcr PT (T) 

dT 



(Ttot exp 



dy 



N 



1 



^-y)Q 2 dk 2 
(i-y) 2 Q 2 



T A q g(a s (k T )) 



Such resummed cross sections improve the perturbative description of differential cross sections, 
dap^/de, for a class of event shapes e, including the thrust. J0| Nevertheless, for small 1 — T, 
that is, close to the limit of massless jets, fits based on Eq. (g) fail to describe the data. This 
is not surprising, because the lighter the jets, the more dependent on long times are their cross 
sections, and correspondingly the more sensitive they are to nonperturbative effects. The mass of 
the jet is a "dial" for tuning the importance of nonperturbative dynamics. 

Although perturbation theory cannot predict true long-time behavior, it can give hints as to 
its nature. In Eq. (|8]) these hints come from the running coupling. When the variable y gets close 
enough to one, the running coupling in the integrals in (|S]) diverges, signalling a breakdown of 
perturbation theory. This divergence is associated with a region of fixed size in kx-, independent 
of N. We should think of this as an ambiguity in perturbation theory, which is resolved in the 
full theory by nonperturbative information. 
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Since the perturbative range of integration 
) should remain meaningful, the minimal modification necessary to (f§) is to replace the 

cto, with 



in Eq. 

lower limit of the fc^-integral with a nonperturbative parameter 



a = 



fi Jo 



dk a s (k) 



(9) 



where fio is a conveniently chosen cutoff. The quantity ao has the interpretation of the integral 
of the running coupling over the nonperturbative region. 

Letting e = 1 — T, we can invert the transform (||[). The substitution (|D then produces a 
simple shift in the perturbative spectrum, 



da(e) da PT (e — X e /Q) 
de de 



+ o 



e 2 Q 2 



(10) 



where we note that the relative size of the effect is X e /eQ, and that corrections begin at l/(eQ) 2 . 
Similar considerations apply to any event shape e that vanishes in the limit of light-like dijets. 

This approach has been formalized and applied in Refs. [b_ 
coupling, Eq. (Of), is thought of as a fundamental, universal parameter. 



The integral of the running 
Applications, which 

include the partial incorporation of higher-order perturbative terms, provide an improved picture 
of differential event shapes, and a unified description of low moments of e = 1 — T. . ., such as 
J dT{l — T)da/dT. In addition, the approximation relates first to second moments: 



X: 



(e 2 ) = (e 2 ) PT + 2-^ (e) PT + ^. 



Q 



Q 2 



This description has been reasonably successful in tying together the first moments of different 
|17| Second moments, however, show large 1/Q 2 corrections, compared to the 

nn 



event shapes. 



predictions of strong-coupling universality, Eq. 
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A more general formalism involves the introduction of nonperturbative "shape" functions |L9 
that generalize Eq. ( |10D to a distribution of shifts due to soft radiation, 



ae Jo ae 



*m ^miz^ + o| _i_ . (12) 



The most important features of this expression are the level of the corrections, down by a full 
power of Q compared to Eq. fllOD , and the independence of the shape function f e from Q. The 
latter implies that a fit to f e at one value of Q, most conveniently Q = mz, is sufficient to predict 
the differential cross section for all Q. fl9| 

Equation (|T2|) may be derived from the factorization properties of soft radiation in perturbation 
theory, and the shape function itself has the interpretation of a matrix element in QCD. To be 
specific, in the case of thrust the matrix element is 

/i-tO, Mir) = (0\W*(p)5 (e - J dn (1 - | cosfl|) £ (n)) W(0)\0) k±<flm , (13) 

with 9 the angle between n and the thrust axis. We define the operators W in terms of path- 
ordered nonabelian phase operators, 



fl^O) = P J9 J °° d\f3-A(\f3) r p e igj~d\(3>.AW) 

and the operators £ are defined to measure energy flow, by pO 



(14) 



N 

£{n)\N) = ^6(cos0-cos0;) 6(<p-(pi) Ei\N) , (15) 
i=i 



for any final state with N particles. The matrix element in Eq. (|T3|) is matched to perturbation 
theory by a cutoff in the transverse momentum of the soft radiation at /xir. In this formulation, 
true universality resides at the level of correlators of energy flow in the presence of the color 
sources, 

Q(n x ...n L]f i) = (01^(0)^(^0 • • • S(H L )W(0)\0) . (16) 

A "mean field approximation" , which eliminates nontrivial correlations between measurements of 
energy flow in different directions: 

L 

Q{ni . . .nx;/i) ->• J\_Q{n h ^i) (17) 

i=i 

reduces the shape function in QTi2| ) to a delta function, / e (e) — > 5(e — A e ). In this approximation, 
the shift ([Top is recovered. 

We may take another viewpoint of Eq. flllf ), and interpret it as a matrix element in an effective 
theory for soft radiation from light-like color sources. This is a natural language for the community 
with special interest in inclusive B decay near the edge of phase space, i.e., with jet-like final 
states, pll] In this case, the color source is the light quark that emerges from the decay b — > sj, 



for example, whose bremsstrahlung factorizes from the remainder of the process. 
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3.3 Exclusive B Decay 



Our third example is the recent application of factorization to the fully exclusive B decays into two 
mesons, Mi, M 2 , either heavy-light {Dir), or light-light (tctt, Kit). Because the incoming quark is 
heavy, the appropriate factorization is somewhat different than in the previous examples, but the 
light-like relative velocity of the light meson(s) in the final state leads to a factorized form for the 
decay amplitude A (here for the light-light case): |^2|, 

A{B^M X M 2 ) = Fb^mJ 1 (m b ,n) ®& Ma (jJ.) 

+f n {m b ,n)®<5> B {fi) ®$ M M ®&Mi(y,fj), (18) 

where the short- distance functions T are computable in perturbation theory, while the $'s are 
nonperturbative matrix elements that are wave functions for the hadrons. Further nonperturbative 
information is contained in Fb-^Mu which is itself a matrix element. It may be possible to compute 
this matrix element if transverse momenta are included in the convolution. [E3| Excitement has 
been generated by the possibility of using the formalism to isolate weak, CP-violating phases in 
these decays. This should be possible because all strong-interaction phases are contained in the 
functions T, to leading power in m&. We still have things to learn about the relationship between 
the different approaches to this factorization, and about the important role of power corrections. 
Nevertheless, the extension of factorization methods to this class of physical problems, whose 
interest transcends QCD, is an important step forward. 



4 BFKL and High Parton Density 

Everything we've discussed so far has involved hard scattering, and hence is restricted to rare 
processes. In recent years, however, considerable attention has returned to the bulk of the high 
energy cross sections, involving relatively low momentum transfers at high energy, including the 
Regge limit (s —>■ oo with t fixed), diffractive scattering and the total cross section. This classic 
constellation of topics is coming to the fore once again, in the light of the copious HERA data on 
small- x DIS, and renewed progress in the perturbative description of the total cross section, the 
so-called "perturbative pomeron" of QCD, as described by the celebrated BFKL equation. 



4.1 BFKL 2000 



The BFKL equation, with the LO kernel shown explicitly, can be written as 

dil)(£,k T ) a s N f d 2 k 



(hr - k' T f 



^,k 7 



V 2 

2 h^f-p 



NLO. 



(19) 



The stimulus for much recent work is the newly-calculated explicit form of the NLO kernel, the 
fruit of a decade of effort. 

As shown in ||25|| , Eq. (fL9|) may be derived from a factorization characteristic of the Regge 
limit. |26| To be specific, we consider the forward scattering amplitude for virtual-photon-proton 
scattering, which is related by the optical theorem to the structure functions of DIS. At low x, we 
need only include (color singlet) gluon exchanges in the t-channel, or equivalently only the gluon 
distribution G(x,Q). Introducing the "unintegrated" gluon distribution ■?/>(£, through 

G(£,Q) = f Q d 2 k T ^,k T ), (20) 
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the relevant factorization for structure function F is 



F(x } Q 2 ) = Jd 2 k T c(^,Q,k T ^^ } k T )(l + 0^ 



x 



(21) 



in terms of a modified coefficient function, c. Relative to Eq. (Q), the roles of kr and £ have been 
exchanged: £ is now the factorization scale, separating, in the terminology of [2l| "fast" from 
"slow" quanta, and kx is now the convolution variable. On the one hand, the incoherence of the 
dynamics of fast and slow quanta make it possible to factorize the amplitude; on the other hand, 
Lorentz invariance leaves the division between the two arbitrary. This arbitrariness leads to an 
evolution equation, the BFKL equation. We note, however, that the factorization in Eq. (f2~T| ) 
holds to next-to-leading logarithm in x ~ s/Q 2 only. Beyond this level, we must generalize the 
equation itself. 



The ansatz ip 



x 



in Eq. flT9|), gives a consistency equation relating the expo- 



nents oj and 7. With a s = N c a s /ir, this is 

^(7) = a s Xo(l) 



1 a a * 1 k ' 
1 -P — ln- 

47r n 



+ "sXi(7) 



(22) 



where the function Xo(j) has long been known, and where xi is new, and the subject of much 
investigation. 

The largest value of to gives the dominant small-x, or equivalently large- s behavior. [24| 
Exhibiting only the LO result in analytic form, one finds 



4AT c ln2(a s /7r)[l - 6.5a s 



tp ~ s 



4iV c ln2 (« s /7r) -large 



(23) 



From LO in the kernel we have QCD Regge behavior, but the innocuous-looking NLO result is, 
as it stands, not quite acceptable. It is simply too large and negative, and can eventually lead, 
not only to a decrease with s, but even to negative cross sections. This produced a bit of initial 
consternation on the part of some enthusiasts, but, ever- resourceful, investigators have developed 
very plausible proposals on how to proceed. In fact, the problem may be traced to "collinear 



:i-7)- 3 . 



which, from the limits k? — > and k? — > k' T in Eq. flTP|), receives poles up 
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divergences' 
to 7~ 3 and 

Proposals on how to interpret the NLO kernel have included: (1) Adjust the scale of a s , 
(2) Impose kinematic constraints in /C, demanding strong ordering of particles in rapidity; |29| (3) 
Import information from DGLAP evolution, given the association of collinear logarithms in DIS 
to the singular behavior of x- [0 Particularly for the latter proposals, the connection of BFKL 
to small-x DIS may suggest phenomenological tests of their efficacy. This story is probably just 
beginning. 



4.2 Effective theories and high parton density 

Sometimes it can be difficult for those not working on small-x and BFKL to appreciate fully their 
perennial fascination. One way of looking at what's special about BFKL evolution is that, if the 
LO BFKL equation is not too misleading, then as we evolve to low x we are forced to a regime 
of high parton density even at "fixed" (actually diffusing) virtuality. One dramatic manner of 
thinking about this regime is as a strong- field configuration of QCD, a dense phase of weakly- 



interacting gluons. pl[ It may even be possible to bring such a state into being in the laboratory; 
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the RHIC at Brookhaven may produce it as an initial state in nuclear collisions, as may the LHC 
operating with nuclear beams. 

This viewpoint has been developed quantitatively through an effective theory, which has some 
similarities to the one described above in the context of shape functions. We introduce a set of 
color sources, this time coming from the distant past along the lightcone, 



W±(x T , x t ) = Pexp 



dx k A*[ 



(24) 



The relation of the BFKL equation to such an effective field theory was described in The 
nuclear connection is made by modelling a large nucleus as a distribution of the sources: f33 | 



^nuclear field = <SqCD + J d 2 X t dx p{x t , X )W+(x , X t ) 



(25) 



Among the intriguing results of this approach is the generation, |34| for a nucleus of essentially 
unlimited size, of a gluon occupation number density, which is of nonperturbative magnitude, 
0(1/ a s ), and which can serve as a starting point for the very complex time evolution of nucleus- 
nucleus collisions. |35] 



5 Conclusions 



Even within the area of factorization at high energy, I have of necessity passed over many de 
velopments from the past few years, regarding global PDF fits 
cross sections at measured transverse momentum, ]38| diffraction 
and skewed 



and their uncertainties, |37| 
9[ higher-twist pEOfl , polarized 
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parton distributions, and more. The subject of power corrections at high 
energy is still new, and we are just now learning to read the quantum mechanical history of QCD 
scattering in the language of final states. I expect the progress of the past few years, punctuated 
as it is with novel ideas and applications, to continue for some time, as we ask new questions of 
quantum chromodynamics. 
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